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Abstract. Using 58 million J/ip decays, we have investigated the pp invariant mass spectrum in the ra- 
diative decay J/%p — > ^ypp and observe a prominent structure with mass near 2m p . Fitting with an S'-wave 
Breit-Wigner, we obtain a peak mass of M = 1859jlf (stat)jl| 5 (sys) MeV/c 2 . J/ip — * 7??c decays from 
the same sample are used to determine the mass, width, and hadronic branching ratios of the r\ c . From a 
sample of 14 million ip(2S) events, the first observation of Xcj (J=0,l,2) decays to AA is made, and branch- 
ing ratios are determined, which are larger than expected from the Color Octet Model. Branching ratios 
of K°K° L in both ip(2S) and J/ip decays are measured, and a more than four sigma deviation from the 
pQCD predicted "12% rule" is observed. In ^(3770) decays, evidence for the non-DD decay to n + ir~ J/ip 
is observed. 

PACS. 13.20.Gd Decays of J/psi, Upsilon, and other quarkonia 
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1 Introduction 

The Beijing Spectrometer (BES) is a general purpose sole- 
noidal detector at the Beijing Electron Positron Collider 
(BEPC). BEPC operates in the center of mass (cm) en- 
ergy range from 2 to 5 GeV with a luminosity at the J/ ip 
energy of approximately 5 x 10 30 cm _2 s _1 . BES (BESI) is 
described in detail in Ref. [1], and the upgraded BES de- 
tector (BESII) is described in Ref. [2] . This paper presents 
some recent results; details may be found in the references. 



The solid curve is the result of a fit using an acceptance- 
weighted S'-wave Breit-Wigner function to represent the 
low-mass enhancement plus the dashed curve to represent 
the background, primarily due to J/ip — > ir°pp where one 
of the photons from the ir° — > 77 is missed. The fit yields 
a peak mass of M = 1859±? (stat)rJ 5 (sys) MeV/c 2 and 
a full width of F < 30 MeV. Here the systematic errors 
include errors determined by generating Monte Carlo sam- 
ples with below threshold peak masses and measuring the 
shift in the output fit masses. 
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2 Studies of J /if) — ► *ypp 

There is evidence for anomalous behavior in the proton- 
antiproton system very near the M p -p = 2m p mass thresh- 
old. The cross section for e + e~ — * hadrons has a narrow 
dip- like structure at a cm energy of y/s ~ 2m p c 2 [3]. In 
addition, the proton's time-like magnetic form-factor, de- 
termined from high statistics measurements of the pp — > 
e + e~ annihilation process, exhibits a very steep fall-off 
just above the pp mass threshold [4]. These data are sug- 
gestive of a narrow, S'-wave triplet pp resonance with J PC = 
1 and mass near M pp ~ 2m p . 

Using 58 million J/ip decays, we have investigated the 
invariant mass spectrum of pp pairs in the radiative pro- 
cess J/tp — > jpp and observe a peak corresponding to 
J/ip — > "fr] c at high mass and a prominent structure with 
mass near 2m p , as shown in Fig. 1. Figure 2(a) shows 
the threshold region for the selected J/tp — * jpp events. 
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Fig. 1. The pp invariant mass distribution for the J/ip 
event sample 
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Fig. 2. (a) The near-threshold M pp 
'ypp event sample. The solid curve is the result of the fit; the 
dashed curve shows the fitted background function, (b) The 
2m p distribution with events weighted by go /<?■ 
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with previous BES measurements and the PDG fit values. 
More detail on this analysis may be found in Ref. [9]. 
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Fig. 3. The m p p invariant mass distribution in the rj c region. 



Further evidence that the peak mass is below the 2m p 
threshold is provided by Fig. 2(b), which shows the M pp — 
2m p distribution when the threshold behavior is removed 
by weighting each event by q /q, where q is the proton 
momentum in the pp restframe and q a is the value for 
Mpp — 2 GeV/c 2 . The sharp and monotonic increase at 
threshold can only occur for an 5*- wave BW function when 
the peak mass is below 2m p . A P-wave BW yields a fit 
which is worse than the S-wetve fit, but still acceptable. 
The structure is not consistent with the properties of any 
known particle. More detail may be found in Ref. [5]. 



3 rj c Parameters 

The mass and width of the r\ c are rather poorly known; 
the confidence level for the Particle Data Group (PDG) 
weighted average mass is only 0.001 [6]. Previously BES 
measured the r\ c mass using the BESI 4.02 million ip(2S) 
sample and obtained M Vc = (2975.8 ± 3.9 ± 1.2) MeV [7]. 
BES also used 7.8 million BESI J /ip events and obtained 
M Vc = (2976.6 ± 2.9 ± 1.3) MeV [8] . For the two data sets 
combined, M Vc = (2976.3 ± 2.3 ± 1.2) MeV and the total 
width r nc = (11.0 ± 8.1 ± 4.1) MeV [8]. 

Here, the mass and width have been determined using 
our BESII 58 million J/ip event sample. We use the chan- 
nels J/ip — > jrjc, with rj c — > pp, K + K~ Tt + n~ , 7r + 7r~7r + 7r~, 
K^KgirT , and (fxf>. Events are selected using particle iden- 
tification and kinematic fitting. Figs. 3 and 4 show the 
mass distributions in the rj c mass region for J 'ftp — > jrj^ 
Tj c — ► pp and rj c — » K + K~ir + ir~ , respectively. Fitting 
the five decay channels simultaneously, we obtain M T/c = 
(2977.5±1.0±1.2)MeVandr^ = (17.0±3.7±7.4) MeV. 
The results for the mass and width are compared with 
previous measurements, including previous BES measure- 
ments, in Figs. 5 and 6. The results are in good agreement 
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Fig. 4. The m K + K - 7T + 7I - invariant mass distribution in the 
T] c region. 

The numbers of r\ c events determined from the fit and 
the corresponding product branching ratios, by decay chan- 
nel, are listed in Table 1. Using the branching fraction 
Br(J/ijj — > 'yr/c) = (1.3±0.4)% [6], preliminary r\ c branch- 
ing fractions can be obtained. Table 2 shows the BES re- 
sults together with the PDG [6] values. The BES Br(r/ C -> 
(f)(f) is smaller than the current PDG value of (7.1 ± 2.8) x 
10~ 3 and is consistent with Belle [10] and DM2 [11] mea- 
surements within errors. Details may be found in Ref. [12]. 



4 X j -> AA 

The lowest Fock state expansion (color singlet mechanism, 
CSM) of charmonium states is insufficient to describe P- 
wave quarkonium decays. Instead, the next higher Fock 
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Fig. 5. Mass measurements of the rj c meson. 
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Table 2. Branching fractions of the rj c . Preliminary. 



Process 


BES(%) 


PDG02(%) [6] 


Br{rj c -> K+K~n + TT-) 


1.2 ±0.4 


2 +U7 


Brijjc — > 7r + 7r _ 7r + 7r - ) 


1.0 ±0.5 


1.2 ±0.4 


Br(rj c -> ^K^) 


1.7 ±0.7 


|(5.5 ± 1.7) 


Br(r)c — » 4>4>) 


0.25 ± 0.09 


0.71 ±0.28 


Br(r] c — > pp) 


0.15 ±0.06 


0.12 ±0.04 



a clear A A signal is observed. After requiring that both the 
7T+J5 and the n~p mass lie within twice the mass resolu- 
tion around the nominal A mass, the AA invariant mass 
distribution shown in Fig. 7 is obtained. There are clear 
Xco, Xci, an d Xc2 — * AA signals. The Jiighest peak around 
the ip(2S) mass is due to ip(2S) — > AA with a fake photon. 
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Fig. 6. Width measurements of the ijc meson. 

Table 1. Number of rj c events and corresponding branching 
ratios for the individual channels. Preliminary. 



Process 

J/4> -* iVc, 


Events 
detected 


Product of 
branching ratios 


rj c — > K + K~ TT + TY~ 


413 ±54 


(1.5 ±0.2 ±0.2) x 10~ 4 


rjc — > 7T + 7r~7r + 7r - 


542 ± 75 


(1.3 ±0.2 ±0.4) x 10~ 4 


rjc -» J R" ± ^7T + 


609 ± 71 


(2.2 ±0.3 ±0.5) x 10" 4 


r; c — > <j5>0 


357 ± 64 


(3.3 ±0.6 ±0.6) x 10~ 5 


r?c -> pp 


213 ±33 


(1.9 ±0.3 ±0.3) x 10~ s 



state (color octet mechanism, COM) plays an important 
role [13, 14]. A calculation of the partial width of Xcj — * PP 
using the COM and a carefully constructed nucleon wave 
function [15], obtains results in reasonable agreement with 
measurements [6] . Generalizing the nucleon wave function 
to other baryons, the partial widths of many other baryon 
anti-baryon pairs can be predicted. Among these predic- 
tions, the partial width of Xcj — > AA is about half of that 
ofxcJ-J>iJ(J=l,2) [15]. 

Using 14 million ipi^S) events and making a scatter 
plot of the 7r + p versus the ir~p invariant mass for events 
with tt + tt~PP mass between 3.38 GeV/c 2 and 3.60 GeV/c 2 , 
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7. Mass distribution of yAA candidates fitted with three 
resolution smeared Breit-Wigner functions and background, as 
described in the text. 



Fig. 



Background from non AA events is estimated from the 
A mass sidebands. The background from channels with 
AA production, including i/)(2S) -> AA, ip(2S) -> S°E°, 
tp{2S) — > AS ± c.c, etc. are simulated by Monte Carlo. 

Fixing the XcO, Xci and Xci mass resolutions at their 
Monte Carlo predicted values, and fixing the widths of the 
three Xcj states to their world average values [6] , the mass 
spectrum (Fig. 7) was fit between 3.22 and 3.64 GeV/c 2 
with three Breit-Wigner functions folded with Gaussian 
resolutions and background, including a linear term rep- 
resenting the non AA background and a component rep- 
resenting the AA background. Fig. 7 shows the fit result. 
The branching ratios of Xcj AA obtained are 
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where the first errors are statistical and the second are 
systematic. 

The results are in contradiction with the expectations 
from Ref. [15], although the errors are large. There is no 
prediction for B(xcO — > AA). More detail may be found in 
Ref. [16]. 



5 Observation of K®K® in if)(2S) decays and 
J /if) decays 

There is a longstanding "p7r puzzle" between J/ip and 
ip(2S) decays in some modes: compared with the corre- 
sponding J/ip decays, many ip(2S) decay channels are sup- 
pressed relative to the pQCD predicted "12% rule" [17]. 
Here ip(2S) — > KgK^ is observed for the first time in the 
BESII 14 million ip(2S) event sample, and the branch- 
ing ratio is used to test the "12% rule" between J/ip and 
ip{2S) decays. 

Candidate events are required to have two charged 
tracks with net charge zero. The two tracks are assumed to 
be 7r + and ir~ , and the fitted intersection of the two tracks 
is taken as the Kg vertex. After requiring the it + it~ in- 
variant mass within twice the mass resolution around the 
nominal K% mass and a Kg decay length in the transverse 
plane longer than 1 cm, the Kg momentum distribution, 
shown in Figure 8, is obtained. Also shown are the Kg 
mass side band events and Monte Carlo simulated back- 
grounds. The signal peak at 1.77 GeV is fitted with a 
Gaussian, and the background below the peak is fitted 
with an exponential function. The preliminary branching 
ratio obtained is (5.24 ± 0.47 ± 0.48) x 10~ 5 . This branch- 
ing ratio, together with results for ip(2S) — > 7r + 7r~ and 
ip(2S) — > K + K~, have been used to extract the relative 
phase between the three-gluon and the one-photon annihi- 
lation amplitudes of ip(2S) decays to pseudoscalar meson 
pairs. It is found that a phase around ±90° can explain 
the result [18]. 

With a similar analysis using the BESII 58 million J / ip 
events, we also measured the branching ratio of J/ip 
KgK^- The Kg momentum distribution in J/ip 
decays is shown in Figure 9. The preliminary branching 
ratio of J/ip -> K S K Q L is (1.82 ± 0.04 ± 0.13) x 10~ 4 . 
This result is significantly larger than the world average 
of ((1.08 ± 0.14) x 10~ 4 ) [6]. Using the branching ratio 
of ip(2S) — ► KgK^, and considering the common errors 
which cancel out in the calculation of the ratio of the two 
branching ratios, one obtains 



K" s Kl 



>h = 



B(i>(2S) 



K° s Kl) 



B{J/if> 



K*Kl) 



= (28.2 ±3.7)%. (1) 



This number deviates from the pQCD predicted "12% 
rule" by more than four sigma. Most interesting is that 
this channel is enhanced in ip{2S) decays, while in almost 
all other channels which deviate from the "12% rule", 
ip(2S) decays are suppressed. These results are prelimi- 
nary. More detail may be found in Refs. [19] and [20]. 




K s Momentum (GeV/c) 

Fig. 8. The K$ momentum distribution in the ip(2S) — » 
K%K\. The dots with error bars are data, the dark shaded 
histogram is from K% mass side band events, and the light 
shaded histogram is the Monte Carlo simulated background. 
The curve shown in the plot is the best fit of the distribution. 




K s Momentum (GeV/c) 
Fig. 9. The K% momentum distribution in J/ip -> K S K° L . 

6 Search for i/?(3770) non-DZ) decay to 

7T + 7T~ J /if} 

The 1^(3770) resonance is believed to be a mixture of 2 3 Si 
and l 3 Z?i states of the cc system [21]. Since its mass 
is above open charm-pair threshold and its width is two 
orders of the magnitude larger than that of the ip(2S), 
it is thought to decay almost entirely to pure DD [22]. 
However, recently some theoretical calculations point out 
that i/)(3770) could decay to non-DZ) final states [23]. 

Here, we report evidence for i/>(3770) — > Tr + ir~J/ip 
based on 8.0 ± 0.5 pb" 1 of data taken in the cm energy re- 
gion around 3.773 GeV with BESII. Another source of 
"7r~ J/ip is from the radiative return process (due to 
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Fig. 10. Scatter plot of the measured 7r + 7r~ energies versus the 
fitted l + l~ masses. There are two clusters. The cluster whose 
energy is around 0.65 GeV is mostly composed of signal events 
from V(3770) — > 7r + 7r~ J/tp, while the events whose energies 
are around 0.57 GeV are due to radiative return events. The 
projections are shown in the middle and bottom plots. The 
open circles are data, the histograms are the results of the 
Monte Carlo simulation for tp(2S) — ► tt + -k~ J/tp, and the solid 
smooth curve is the fit to the data. 

initial state radiation (ISR)) to the tp(2S) followed by 
tp(2S) — > ir + ir~J/ip. We developed a new generator is- 
rpsi which includes production of J/tp, tp(2S) and other 
resonances due to radiative return. The Monte Carlo sim- 
ulation includes leading-log-order radiative return, where 
the cm energies after ISR are generated according to Ref. 
[24]. 

To search for the decay of ^(3770) — > 7T+7T J/tp and 



J/tp 



or [i + n 



and [V \x 



can- 



didate events are selected. They are required to have four 
charged tracks with zero total charge. 

Figure 10 shows a scatter plot of the 7r + 7r~ energies 
versus the fitted masses of the l + l~ after a four constraint 
fit to the process ^(3770) — ► 7r + 7r~ J/tp. In the middle 
sub-figure, the higher mass peak is due to radiative return 
to the ip(2S) followed by tp(2S) -> nirJ/tp. This peak is 
shifted to 3.18 GeV because E cm is set to 3.773GeV in 
the kinematic fitting. The fit to this peak yields a total 
of 2.2 ± 0.4 background events near the signal peak at 
3.097 GeV, out of the 9.0 ± 3.0 events. After background 
subtraction, 6.8 ±3.0 signal events remain. The branching 
fraction for the non-DD decay V(3770) — > n + n~J/tp is 
measured to be 



BF(tP(3770) 



-J/tp) = (0.59 ±0.26 ±0.16)%, (2) 



where the first error is statistical and the second system- 
atic. Using the total width of the ^(3770) resonance from 



the PDG [6], this branching ratio corresponds to a partial 
width of 

1^(3770) -> 7T+7T- J/tp) = (139 ± 61 ± 41) keV. (3) 

More detail may be found in Ref. [25]. These results are 
preliminary. 



References 

1. J. Z. Bai et al, (BES Collab.), Nuc. Inst. Meth. A344, 319 
(1994). 

2. J. Z. Bai et al, (BES Collab.), Nuc. Inst. Meth. A458, 627 
(2001). 

3. A. Antonelli et al, Nucl. Phys. B517, 3 (1998). 

4. G. Bardin et al, Nucl. Phys. B411, 3 (1994). 

5. J. Z. Bai et al, (BES Collab.), Phys. Rev. Lett. 91, 022001 
(2003). 

6. K. Hagiwara et al, Phys. Rev. D66, 010001 (2002). 

7. J. Z. Bai et al, (BES Collab.), Phys. Rev. D60, 72001 

(1999) . 

8. J.Z. Bai et al, (BES Collab.), Phys. Rev. D 62, 72001 

(2000) . 

9. J.Z. Bai et al, (BES Collab.), Phys. Lett. B555, 174 (2003). 

10. H.-C Huang et al (Belle Collaboration), hep-ex/0305068. 

11. D. Bisello et al. (DM2 Collaboration), Nucl. Phys. B350, 
1 (1991). 

12. J. Z. Bai et al, (BES Collab.), submitted to Phys. Lett. 
B, hep-ex/0308073 (2003). 

13. See, for example G.T. Bodwin, E. Braaten and G.P. Lep- 
age, Phys. Rev. D51, 1125 (1995); Han-Wen Huang and 
Kuang-Ta Chao, Phys. Rev. D54, 6850 (1996); J. Bolz, 
P. Kroll and G. A. Schuler, Phys. Lett. B392, 198 (1997). 

14. J. Z. Bai et al. (BES Collab.), Phys. Rev. Lett. 81, 3091 
(1998). 

15. S. M. Wong, Eur. Phys. J. C14, 643 (2000). 

16. J.Z. Bai et al, (BES Collab.), Phys. Rev. D67, 112001 
(2003). 

17. J. Z. Bai, et al. (BES Collab.) Phys. Rev. D67, 052002 
(2003). 

18. C.Z. Yuan, P. Wang and X.H. Mo, Phys. Lett. B567, 73 
(2003). 

19. J. Z. Bai et al, (BES Collab.), submitted to Phys. Rev. 
Lett., hep-ex/0310024 (2003). 

20. J. Z. Bai et al, (BES Collab.), submitted to Phys. Rev. 
D, hep-ex/0310023 (2003). 

21. P.A. Rapidis et o/.(MARK-I Collaboration), Phys. Rev. 
Lett. 39, (1978) 526. 

22. W. Bacino et al, Phys. Rev. Lett. 40, 671 (1978). 

23. H.J. Lipkin, Phys. Lett. B179, (1986) 278. 
Y.P. Kuang, Phys. Rev. D65 094024 (2002). 

24. E.A. Kuraev and V.S. Fadin, Yad Fiz, 41, 377 (1985). 
(Sov. J. Nucl. Phys. 41, 466 (1985)) 

25. J.Z. Bai et al, (BES Collab.), hep-ex/0307028 (2003). 



